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Abstract 

Since the discovery of the large angular scale anisotropies in the microwave background radi- 
ation, the behaviour of cosmological perturbations (especially, density perturbations and gravi- 
tational waves) has been of great interest. In this study, after a detailed and rigorous treatment 
of the behaviour of gravitational waves in viscous cosmic media, we conclude that the damping 
of cosmological gravitational waves of long wavelengths is negligible for most cases of physical 
interest. A preliminary analysis suggests that similar results hold for density perturbations in 
the long wavelength limit. Therefore, long wavelength cosmological perturbations have not been 
practically affected by viscous processes, and are good probes of the very early Universe. 
Key Words: Shear Viscosity, Damping, Viscous Cosmic Media, Amplitude of Gravitational Waves, 
Density Perturbations 

1 Introduction 

A likely explanation of the observed large angular scale anisotropy of the Cosmologigal Microwave 
Background Radiation (CMBR) is cosmological perturbations of quantum mechanical origin (Gr- 
ishchuk 1993), mainly gravitational waves (Grishchuk 1994). However, one could argue that these 
perturbations are not responsible for the CMBR anisotropy if they could be washed out by viscous 
processes. As a first step, one should investigate whether, in a classical approach, the amplitude of 
such perturbations decreases significantly or not. This is the point of this study. We have dealt mainly 
with gravitational waves in several viscous cosmic media, posscsing only shear viscosity. One could 
refer to some standard sources, like Hawking 1966, Weinberg 1972 and Grishchuk and Polnarev 1980. 
We adopted a technique previously used by Weinberg 1972. Several particular models for viscous 
cosmic media were taken from Mendez et al 1997. 

The assumptions used in the study of the early Universe can be found in standard text books, 
for example Weinberg 1972. Our study indicated that the zero chemical potential approximation 
should be reconsidered by the end of the radiation dominated era, in order to get correct decoupling 
time between matter and radiation, and the hydrogen and baryon abundances should be taken into 
account. We neglect the bulk viscosity. 

Information concerning transport phenomena can be found in Landau and Lifshitz 1966, Tabor 
1970 and Reif 1965. We will only include the information necessary for our study. By "viscosity" we 
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t^Th 


Perfect thermal equilibrium holds. Dissipation is negligible. 


t<Th 


Departures from thermal equilibrium begin. At the equality sign, decoupling occurs and 
this interaction is no longer realised. Dissipation becomes important ,untill of course 

decoupling occurs and the interaction ceases. 


t>Th 


Decoupling has occured, and no dissipation takes place since the interaction has stopped. 



Table 1: LOCAL THERMAL EQUILIBRIUM, DISSIPATION AND DECOUPLING 



will mean specifically shear viscosity. If n is the mean free time between collisions, et.i is the thermal 
energy density of the particles of the i species responsible for momentum transfers, then the shear 
viscosity coefficient is 

^ = OY^et.in (1) 

i 

where ^ is a numerical constant. 

Elementary information required about Local Thermodynamical Equilibrium, Dissipation and De- 
coupling are presented in Table 1 (Weinberg 1972, Kolb and Turner 1990). We should mention though 
that these conditions are a consequence of the constancy of the entropy per comoving volume element. 
We define a characteristic Hubble time by Th = TZ/TZ where TZ is the scale factor, t] is conformal 
time, a dot is differentiation with respect to t, while a prime with respect to rj. 

2 The behaviour of Gravitational Waves in a non viscous 
medium 

Lifshitz pointed out how the different types of perturbations can be constructed in the form of scalar, 
vector and tensor harmonics, corresponding to density perturbations, rotational perturbations and 
gravitational waves respectively (see, for example, Lifshitz and Khalatnikov 1963). The perturbed 
Einstein equations are given in Weinberg 1972. For gravitational waves, hj denotes the time dependent 
part of the perturbations and q the constant wavenumber of the perturbation, related to the time 
dependent wavelength A(t) by the relation q = Since all the h-^ components obey the same 

equations, we will ignore the indices in our notation and refer to a single component, as h. This 
equation, has been transformed and interpreted as a parametrically excited oscillator. (Grishchuk 
1993). The /i(?7) amplitude is related to the h(rj) amplitude by ^,(77) = ii{ri)Ti{ri)^^ . The TZ~^ 
variation of h reflects the adiabatic decrease of h. The perturbations interact with the background 
time dependent gravitational field, which supplies energy to waves with wavelengths that satisfy 
the parametric amplification condition. In the case of gravitational waves, the interaction potential 
U{rf) = represents the background gravitational field. In order for this interaction to take place, 

the frequency of the wave must be comparable with that of the variations of the background field. 
Depending on the wavelengths of the perturbations, the behaviour of fj, and h is as described in Table 
2. Therefore, gravitational waves interact parametricaly with the background gravitational field, and 
the ones longer than the Hubble radius are "superadiabatically" amplified. The quantum treatment 
of this phenomenon and its implications to the CMBR statistics and anisotropy is investigated in 
Grishchuk 1993. 
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In this overbarrier region, the osciUatory solutions 
experience only adiabatic decrease, since no 

1 nfprP pt.imi witn tnp njirripr tPiVp^ nl^tpp T'np 

expressions for the amplitudes are 
^(77) = Ci exp{±iqri} and 
h{r]) = CiTZ"^ exp{±j(777}. 


g«:C/(?7) 


In this underbarrier region, the solutions are practically 
constant. These long wavelengths do not suffer from adiabatic 
decrease: due to their interaction with the barrier they are 
amplified (" superadiabatic amplification" ). The expressions 
for the amplitudes are 
Ai(ry) = CbTl + Call J n-^d{r]) ~ CbTZ 
and h{7j) ^Cb + CaJ TZ-^dirj) ~ Cfc. 



Table 2: AMPLITUDE OF GRAVITATIONAL WAVES IN THE NONVISCOUS CASE 



3 The behaviour of Gravitational Waves in the presence of a 
Shear Viscosity 

In order to put the equations of propagation for gravitational waves in a viscous medium in a form 
that reveals the underlying physics, we will first derive an expression for shear viscosity. We begin 
from equation (|l|). It is obvious that species with Ct.iTi several orders of magnitudes smaller than those 
of the rest constituents of the fluid will not participate. We will make some assumptions. The first is 
that in the radiation dominated era, the thermal energy density of a species is approximately equal 
to its total energy density: et.i — ti. The second assumption is that the ti of the different species are 
of the same order of magnitude. We define the number ki = e^/e, where epsilon is the total energy 
density of the fluid. The different will be of the same order of magnitude, and will all be denoted 
by K. Therefore, (0) reduces to 

i = ei^n,Ti)e (2) 

i 

Obviously, the quantity f = t^iTi) is the mean free time of the fluid, and will be at most of order 
of the largest ti. The third assumption, is that the Ti of the species contributing to ^ are of the same 
order of magnitude and equal to t. Then, (^) becomes 

C = e^Ker (3) 

The fourth assumption is that the number of species contributing to viscosity is small: ■0 1- 

Following Weinberg 1972, one finds that in the presence of a cosmic medium possesing shear 
viscosity the propagation equations for the gravitational waves' amplitude is 

/i+(3T^i+66lV'«:T^V)/i+(^)2/i = (4) 

where we have used and expressed e in terms of from the unperturbed Einstein equations. 
Equation (^) shows that a further (above adiabatic damping) decrease of the amplitude arises due to 
shear viscosity. This equation in terms of conformal time and for the /i amplitude is 

/x" + 66lV'KA2^/i' + (g2_y(r;))^ = (5) 
where A = ^ and the time dependent potential is V{r]) — [/(ry) + GOi/jkA^^. We put (||) into the 
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form of the Schroedinger equation by introducing the function m: 

fj, = mexp{- J SeijuA^^difj)} (6) 

Then, equation ^ becomes 

m" + (q^ - Y{r]))m = (7) 

with Y{ri) = F' + F'^ where F ^ A + ^OipnA^^. This is again the equation for a parametrically 
excited oscillator, but with a modified potential due to shear viscosity. 

It is known that the dissipation is expected to be negligible when the mean free time between 
collisions is much less than the Hubble time. This is a consequence of the constancy of entropy within 
a comoving volume element, which yields the relation TZ oc T~^. The expansion rate T^^ determines 
the rate of temperature change: particles that are in thermal equilibrium should have an interaction 
rate greater than the rate of temperature change. For r oc T^^, like the ones we will consider, when 
T > Tff, a particle will interact less than once, so this species will drop out of equilibrium. 

When r < Th, the viscous term participating in the cofactor of /i in (|4|) is much less than the 
expansion term. Even in the case of non negligible viscosity, the viscous term can never become larger 
than the expansion term, because otherwise the condition of thermal equilibrium (not necessarily 
perfect thermal equilibrium) (see Table 1 ) will be violated. This means that decoupling will occur and 
the dissipative mechanism under consideration will cease to function as such. Significant dissipation is 
expected to occur at those times when the viscous term becomes comparable to the expansion term: 
then, we have large departures from perfect thermal equilibrium. 

Let h he a. solution in the absence of viscosity and h a solution of (^, that is, in the presence of 
viscosity. A measure of the dissipation can be presented as 

h~h 

Z = (8) 
h 

We will now consider short and long wavelengths separately. 
3.1 Solution for short wavelengths 

A solution of (^ in this limiting case of short waves is (Weinberg 1972) 



h = /iexp{- j WipnT^^rdt} (9) 

This solution is applicable under the condition Th ^ that is, when the relevant waves are well 
inside the (time dependent) Hubble radius. This result can be obtained from (^ by neglecting the 
potential Y[ri). Combining (^ with (^, one derives for this limiting case the damping 

Z =1- cxp{- J WipnT-'^Tdt} (10) 

There is no a priori reason for this damping to be much smaller than 1, nevertheless, it will be at 
most of order unity, otherwise, the condition of thermal equilibrium will be violated. 

3.2 Solution for long wavelengths 

This is the limit which, apparently, was not considered before for viscous matter. In this case, 
wavelengths are much longer than the Hubble radius in the era under consideration and, therefore, 
the term in (Q) can be neglected. We have found a solution of equations (Q) and (Q) in the form 

h^Ci+C2 j TZ-^{exp{-2 J Wxl}nT^^Tdt})dt (11) 
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The solution in the nonviscous case is h — Ci + C2 J TZ~^dt so viscosity affects only the "decaying" 
second term, but not the constant ("growing") first term. Even in the absence of viscosity, the 
"decaying" term is much smaller than the "growing" one. Thus, since the affected by viscosity 
"decaying" term in h is smaller, or, at most, equal to the "decaying" term in h, it will also be much 
smaller than the (unaffected) constant term. Therefore, 

This means that the viscosity does not practically affect the amplitude of waves with wavelengths 
longer than the Hubble radius. Of course, the same holds for the wavelengths that are longer than 
the Hubble radius today. 

This is a central result of our study, because these long wavelengths are responsible for the large 
angular scale CMBR anisotropy. In this case, we derive for the absorption the expression 

C2(/ n-^dt - / 7^-3(exp{-2 / Se'iPK,Tjj^Tdt})dt) 
As we have already argued, this is always much smaller than 1. 



4 Damping of short waves in Various Viscous Cosmic Fluids 

To derive concrete numbers, our strategy is to compare r with Th and then calculate Z defined by 
([To|). Th was taken as Th — 2t, and r, Th, and Z were expressed as functions of the temperature T. 
We are dealing with radiative fluids, where viscosity arises due to failure of perfect thermal equilibrium 
between matter and highly relativistic particles, like photons and neutrinos. The mean free time was 
calculated as the inverse of the product of the velocity of light with the relevant cross section, a, and 
with the particle density n of the particles of matter which interact with the radiation: t = {anc)~^ . 
In the case of Quark Gluon plasma, the expression of the mean free time was taken from Thoma 
1991. All the scenarios are realised in the Radiation Dominated Era. tpi denotes the Planck time, 
"nipi, mp, me, m^^ the Planck, proton, electron, muon masses respectively, a the fine structure constant, 
k the Boltzmann constant while L is the product of the reduced Hubble parameter with the baryon 
density parameter: L ~ flBh^- is a dimensionless number between 0.4 and 1 that represents the 
uncertainty to the observed value of Tj^^ due to systematic errors, fi^ is the ratio of the baryon 
density of the Universe over the critical density. We have adopted for L the value 2.5 • 10~^. 
The time — temperature relation is given by 

t = 0.3(50-^/^(^)%. (13) 

(see Kolb and Turner 1990) where 5, denotes the relativistic degrees of freedom (number of effectively 
massless degrees of freedom, mc^ ^ kT), and varies with time. 



4.1 Quark Gluon Plasma 

According to Mendez et al 1997, the period of interest is (10^"^ > T > W^)K. Taking the expression 
of the mean free path from Thoma 1991, the mean free time of quarks is given by 

r, = 3.8.10^V (14) 

reducing to the expression 

T, = 1.44- 10-20 ilL£s (15) 
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Similarly, for gluons 
reducing to the expression 



Tg = 1.3 ■ W^^^tpi (16) 



T 

The relation between time and temperature reduces to 



r<, = 4.83 -10-21^^5 (17) 



t = 3.2.10-\^rs (18) 

5* is 106.75 for the quark giuon plasma era. By considering 3 generations of quarks and 8 kinds 
of gluons, we found k = 10^^, ^ — 14 and used r 10~^° ^^^^ s. We took 6 — 4/15 after 
consulting Thoma 1991. Then, Z - 4.13 ■ 10"^ Equation (|) gives a shear viscosity coefficient 
^ — 7 ■ 10~^^K~^T^ g cm^^ s^^, while the viscosity coefficient resulting by a straightforward use of the 
expression of ^ found in Thoma 1991 is ^ = 4 • \Q~^'^K~^T^ g cm^^ s^^ . There is a good agreement 
(the ratio of the shear viscosity coefficient of Thoma 1991 to our shear viscosity coefficient is at most 
5). The absorption due to this mechanism is small, but there is a significant departure of our result 
from the one given in Mendez et al 1997, which is Z ~ lO""^. This can be explained from the fact that 
in Mendez at al 1997 the shear viscosity coefficient used is given as ^ = 1.88 • IQ^^^ K^T^ g cm~^ s~^. 
For the derivation of the latter, Mendez et al 1997 also quote Thoma 1991, and they use the same 
assumptions as we do. 

An interesting feature of this absorption mechanism is that the quantity AT^^t (thus, the shear 
viscosity coefficient) is proportional to T^, decreasing with time! This is the only considered medium 
having this behaviour. Large deviations from perfect thermal equilibrium occur in the beginning of 
this era and not its end, unlike the other considered media. Those deviations decrease and the fluid 
approaches perfect thermal equilibrium in the course of time. Viscosity is important towards the 
beginning of this era, when the density is large and the mean free time is small. As can be seen from 
(|l|), the density and the mean free time are the parameters to determine viscosity. In the course of 
time, the decrease of density is accompanied by increase of the mean free time. Depending on the fluid, 
either one or the other primarily determine the viscosity. Since viscosity in the quark gluon plasma 
is greater at large densities, one expects that the collisions that cause the transport of momentum 
do not involve great momentum transfer between particles, but frequently take place and cause the 
dissipation to occur. 



4.2 Electron Neutrino Mixture 

Two kinds of neutrinos are of main concern for this model, the muon and the electron types, and 
their antineutrinos. Muon type neutrinos decoupled when the muons annihilated at T ~ 1.2 • lO^^K, 
since their reaction rate is sensitive to the presence of muons (Weinberg 1972). Electron type neu- 
trinos decoupled later, when T ~ (time of electron — positron annihilation) because their re- 
action rate is sensible to the presence of electrons. Concerning electron and muon type neutrinos 
and antineutrinos, electrons, positrons and photons to contribute in the energy density of this fluid, 
Ky = = = Ky = Kcr = 10~^. = 4/15, in agreement to Weinberg 1971. 

The period of interest is (10^^ > T > 10^'^)K. Taking into account the cross section of weak 
interactions and the electron density (de Groot ct al 1980), the mean free time of electron type 
neutrinos and antineutrinos is given by 

r.,=2.8.10"(^)4(^)%, (19) 
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reducing to 
For muon type neutrinos and antincutrinos 



T,, = (1.47-10^^)^s (20) 



r., = 3.9 . lO'"i^)H^)' eM^}t,i (21) 

reducing to 

1 O^O 1 9'^ • 1 K 

r.,. = 2.1 . 1010(1^)5 exp{l:^iJ£-^} (22) 

The time — temperature relation is 

t = l.lCy^fs (23) 

The laXio t^Jtu, = 1.41 ■ exp{1.23 ■ lO^^K/T} increases rapidly in the period (lO^^ >T> 1.2 -10")^ 
in a range of half up to three orders of magnitude. Therefore, we shall study the behaviour of the 
fluid in two subperiods: 

4.2.1 Subperiod (lO^^ > T > 1.2 • 10")ii: 

Due to the rapid increase of Ti, / t^^ , it is more convenient to express the shear viscosity coefficient 
through (^ and compare the average mean free time 

2 

in c 

f = 2K,r,^(l + 1.41exp{-^}) (24) 

with Th- The contribution of muon type neutrinos and antineutrinos to viscosity is greater than that 
of their electron type counterparts. This is due to the faster growth of t^^ in this subperiod. The 
shear viscosity coefficient results to be 3.17 • W'>K{1 + 1.41 exp{1.23 • W^K /T})T-^g cm-'^ s^^. 



4.2.2 Subperiod (1.2 • lO" > T > IQ^°)K 

From now on, only electron neutrinos and antineutrinos contribute to the shear viscosity. We use (H) 
with 1^ = 2 and the mean free time t^^. The resulting shear viscosity is 3.17- K g cm~^ s~^, in 
good agreement to the shear viscosity resulting from deGroot et al 1980, which is 2.68-10^^7^ T'^^g cm~^ s 

Dissipation is Z ^ 4.11 • 10~^ from both periods. The muon neutrino and antineutrino contribution 
turns to be of the same order, but slightly larger than that of their electron type counterparts. This 
is due to the faster growth of r^^ being compensated by the shorter time period of contribution of 
the muon type neutrinos and antineutrinos to viscosity. We are not aware of anyone having reached 
a similar or a contradicting result. Neglection of the muon neutrino and antineutrino contribution, 
leads to Z^^ - 6.19 • IQ-^. 

Mendez et al 1997 use for the shear viscosity coefficient the expression f — 8.79-10^'^i4rT~^ g cm^^ s^^ 
(they refer to deGroot et al 1980 as the source of their derived shear viscosity coefficient as well). This 
can explain a departure of two orders of magnitude between their result and ours, but cannot justify 
the actual departure of seven orders of mag nitude. Indeed, Z 5 • 10"^ in Mendez et al 1997. These 
authors have not considered the contribution of the muon neutrinos, but this can not explain the 
difference. 
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4.3 Thomson Scattering 

The period of interest is (10^ > ^ > 3 • 1Q'^)K (Padmanabham 1993). The mean free time in Mendez 
et al 1997 is a result of the zero chemical potential approximation (also applied to the previous 
mechanisms). Nevertheless, in these stages of the radiation dominated era, prior to recombination, 
this assumption seems to be not valid, and the hydrogen and baryon abundances should be taken into 
account. Indeed, setting this mean free time equal to the Hubble time, one will reach the result that 
decoupling occured at some temperature between 10^ K and 10^ A'. Following up to a point Kolb and 
Turner 1990 (the Thomson cross section and number density were taken from Weinberg 1972), we 
adopted the mean free time 

' (25) 



-l + (l + C2^3"^/'e'5^)l/2 



With Ci = 7.56a-2(f-)i/2, C2 = 10"^ S, = ^2 = ""''' ^ ^3 = (H reduces to 



"■nipi' ' ^" J '-'I kT ' kT ' "■i kT 

4.2 • 10-2j^3/2(r)-3/2g^p( 1.6-10^g ) 

-1 + (1 + 2.2 ■ 10-22X-3/2(r)3/2exp( i-6-^o"-^ ))"V^ 



(26) 



For (3.4 • 10^ > T > 3 • 10^) AT, our mean free time reduces to r = 2.8 • lO^Ar^Z-^T-^/^ exp{7.9 • 
10^ A:T~1}s, in good agreement to t = 8.7 • lO^AT^/^T-^/'* exp{8 • WKT-^}s which holds for tem- 
peratures close to photon decoupling (Padmanabham 1993). 

The time — temperature relation is reducing to 

t^2C^fs (27) 

Absorption is more significant towards the end of the period of consideration [T < 3.6 • lO^if), 
when departures from perfect thermal equilibrium become large. The main contribution is from this 
subperiod. k = 7 ■ 10~^, ip = 1 and = 4/15. We have calculated Z ~ 7.14 • 10^"'. Our result agrees 
with that of Mendez et al 1997 , who give Z ^ 2 ■ We cannot explain this agreement in results, 

since we have adopted completely different assumptions. We abandon the zero chemical potential 
approximation, and in Mendez et al 1997 this mechanism is considered only up to 4 • lif'K . Their 
adoption of zero chemical approximation gives physically wrong results (decoupling of photons between 
lO^if and l{fK, and also a damping due to Thomson scattering six orders of magnitude greater than 
for Compton scattering, even though the latter is a more efficient thermalising mechanism). Besides, 
their shear viscosity coefficient 8 • lO^^^AT^^/^T^/^ exp{6.02 • lO^ATT"^} could never lead to a result 
such as theirs, because the exponential blows up towards the end of the era of applicability of Thomson 
scattering, even within the temperature range adopted in Mendez et al 1997. 

4.4 Compton scattering 

The applicability of this mechanism is within (10^ > T > 5.8 • 10'^)K. Using the cross section given 
in Padmanabham, the mean free time turns out to be 

' -t,i (28) 



-l + (l + C2^3"^/'e'52)l/2 

where C3 = 7.56q!~^(^^)'^/^. This expression is greatly simplified to 

r « 2 • lQ^°K'^T-^s (29) 
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after a Taylor expansion of the square root in the denominator. This is in good agreement to r = 
2.17 • lQ^^K'^T~^s (Padmanabham 1997). Our shear viscosity coefficient is ^ = 5.6 • 10^'^ g cm~^ s~^, 
roughly two orders of magnitude smaller than the ^ = 4.37 • 10^^ g cm~^ derived by using the 
mean free time of Padmanabham 1993. The result is Z 6 • 10~^, presenting a departure of seven 
orders of magnitude from the Z ~ 10~^ of Mendez et al 1997. This can be explained by their use of 
^ = 1.72 ■ 10^ g cm-'^ 5"^ (they quote Padmanabham 1993). 

The damping due to this mechanism is two orders of magnitude greater than that of Thomson 
scattering. This is in agreement with the fact that Compton scattering is a much more efficient 
thermalising mechanism than Thomson scattering (Padmanabham 1993). 

5 Damping of Density Perturbations 

The treatment of density perturbations in the non viscous case (for example Grishchuk 1994) shows 
that they interact with the background gravitational field in the same manner as gravitational waves. 
Since one expects that viscosity will always act against the deformation of the medium produced by a 
perturbation independently of the nature of the perturbation, one expects that density perturbations 
should have a behaviour similar to gravitational waves in presence of viscosity as well. In fact, this 
should be true for a certain range of wavelengths, since we know that density perturbations of smaller 
wavelengths are washed away. But in the longer wavelength limit one could still expect that the 
behaviour of density perturbations and gravitational waves should be quite similar. Of course, this 
issue requires a more rigorous treatment, following the same steps we have already performed for 
gravitational waves. 

6 Discussion of results 

The main result of this study has been that gravitational waves which are longer than the Hubble 
radius today are not practically affected by dissipation, and, consequently, there has been no damping 
in their amplitude. Thus, these modes are good candidates for the production of the CMBR large 
angular scale anisotropy, although a study at the quantum level is required in order to investigate 
whether viscosity affects more delicate properties such as squeezing. 

For shorter waves, the picture is different. These modes have been affected by damping more 
severely than the longer ones. The absorptions we calculated for shorter wavelengths are given in the 
table of results. Compton scattering is the most efficient damping mechanism, Thomson scattering 
and viscosity due to electron neutrino mixture are of comparable efficiency, while the effect of quark 
gluon plasma in the damping of gravitational waves is much smaller. In the electron neutrino mixture, 
it seems that the muon type neutrinos and antineutrinos play a slightly more important role than the 
electron type ones (they both give same order of magnitude damping, with the numerical coefficients 
slightly in favour of muon type neutrinos and antineutrinos). 

Thus, the more time a wave has been within the Hubble radius, the more damped will its amplitude 
be. The ones longer than today's Hubble radius are not affected at all. 

There is a large departure of our results from those of Mendez et al 1997, possibly due to numerical 
disagreements. Our estimates of damping are several orders of magnitude higher than the ones of 
Mendez et al. According to these authors, these results are upper limits of damping, and are given in 
Table 4. In the case of Thomson Scattering, they use the zero chemical potential approximation, which 
as we have shown should be reconsidered, and apply this mechanism for a shorter period of time. In 
the case of other mechanisms, the source of our disagreement lies elsewhere. We have used the same 
references for the expressions of mean free paths, cross sections and number densities, and the same 
assumptions. Their approach is much more complicated than ours: they set the Einstein equations in 
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Viscous Cosmic Medium 


Era of Application 


Damping 


Quark Gluon Plasma 


(10^^ >T> 10^"*)/^ 


Z ~ 1.55-10-'' 


Electron Neutrino Mixture 


(10^^ > T > 10^")i^ 


Z ~4.11 • 10"^ 


Thomson Sc;attcring 


(10" > T > 3.1 • W)K 


Z - 7.14-10-'' 


Compton Scattering 


{10^ > T > 5.8 • 10'^)K 


Z~ 6.14 -10-^ 


Table 3: TABLE OF RESULTS 




Viscous Cosmic Medium 


Era of Application 


Damping 




Quark Gluon Plasma 


(10^^ >T> 1Q'-^'^)K 


^10-^ 


Electron Neutrino Mixture 


(10^^ >T> 10''>)K 


^5-10-y 


Thomson Scattering 


(10'-* > T > 4 • 10''^)K 


- 2 • 10-^ 


Compton Scattering 


(10" > T > 5.8 • lO'');^' 


~io-y 



Table 4: TABLE OF RESULTS OF MENDEZ ET AL 1997 

quasi-Maxwellian form and retain in their formulas the shear viscosity coefficient. Then, they solve the 
resulting differential equations either analytically or numerically. If one calculates the shear viscosity 
coefficients directly from the references, these turn out to be several orders of magnitude larger than 
those used in Mendez et al 1997. These departures can explain the differences of our results with the 
results of Mendez et al 1997. 
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